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Abstract 

The gap function <ft(k,u), determined from a Lanczos calculation for a doped 2-leg t-J ladder, is 
used to provide insight into the spatial and temporal structure of the pairing interaction. It implies 
that this interaction is a local near-neighbor coupling which is retarded. The onset frequency of the 
interaction is set by the energy of an S = 1 magnon-hole-pair and it is spread out over a frequency 
region of order the bandwith. 
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The momentum and frequency dependence of the superconducting gap function (f)(k, uj) 
provides information on the spatial and dynamic structure of the pairing interaction. Tra- 
ditionally, electron tunneling has been used to explore the frequency dependence of the 

n n 

gap Presently, angular resolved photoemission spectroscopy |2j (ARPES) and scanning 
tunneling microscopy |3j, |4| (STM) open the possibility of obtaining both k and uj informa- 
tion about the superconducting gap. However, approaches to extract this information on 
4>(k, uj) from such experiments are still being explored. As has been previously discussed, a 
k and uj dependent gap function for finite t-J lattices can be obtained using Lanczos exact 
diagonalization j^, . Here, using results for <p(k, uj) obtained for a doped 2-leg t-J ladder, 
we explore what can be learned once data for (p(k, uj) becomes available. 
The Hamiltonian for a 2-leg t-J ladder can be written as 

H = J rung ^ ] ( Sn ■ S{2 — — nnrii2 J + Jleg ^ ] f Sj + i iQ ■ Sj jQ — — J 

i ^ ' i,a ^ ' 

- Wg ( C ll C ^2 + h.C^j - t lcg (4+l,a C ia + h.C^j (1) 

i i,a 

where Q a are projected hole operators (spin indices are omitted) and a = (1, 2) labels the 
two legs of the ladder. We will consider the isotropic case in which J rung = J\ cg = J and 
our calculations will be carried out for a periodic 2 x 12 ladder at 1/8 
and 1/6 hole doping. 

As previously discussed, a gap function can be extracted by combining Lanczos results 
for the usual one-electron Green's function G(k,u) with the Fourier transform of Gorkov's 
off-diagonal Green's function 

F(k, t) = i (Tc_ fcl _ ff (t/2)cv(-*/2)) (2) 

Here, for a finite system, this expectation value is taken between the ground states for N and 
N — 2 particles and we choose the phase of F to be zero. We will take iV = 22, iV — 2 = 20 
corresponding to an average filling n ~ || = 0.875. The Dyson equations j?J relating G and 
F are 

[Z{k, ui)u- {e k + X(k, to)] G(k, to) = 1 - <f)(k, uj) F(k, uj) (3) 
[Z(-k, -uj) uj-(e k + X(-k, -uj))} F(k, uj) = uj) G(k, uj) (4) 
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with Z and X the usual Nambu self energies and <p(k, u) the gap function. Then, using the 
even {k — > —A;) parity of Z(k,u), X(k,u),€k and <p(k,uj) along with the even u dependence 
of Z(k,u), X(k,u>) and (p(k,u), we have 

\ F(k,u) 

W' u >- F*(k,u>) + G(k,u)G(k,-u)- [ ) 

Alternatively, Z(k,u) can be eliminated to obtain an expression for the superconducting 

gap 

l ' } Z(k,co) G{k,u)-G(k,-u)' { ) 

From a numerical calculation of G(k,u) and F(k,u), A(k,u) has been obtained for a 32- 

site, t-J cluster jjj. This extended previous work by Ohta et. al who fit the spectral 

weight Im F(k,u) to a d x 2_ y 2 BCS-Bogoliubov quasiparticle form in which the frequency 

dependence of the gap was neglected. 

In contrast to the long range order of the superconducting ground state of a 2D lattice, a 

2-leg ladder exhibits power law pair field correlations 11] which decay as x~ l l Kp . Here, 

k is the Luttinger liquid parameter associated with the massless charge mode. This implies 

nnn 

that for a ladder of length L, the off-diagonal Green's function F(k,ou) decays [12J, [13J, Il4j 

i 

as (£ / L) 2k p . Here, the coherence length £ is proportional to the inverse of the gap. For our 
doped ladder L = 12 is of order this coherence length so that we can probe the internal 
k — to structure of a pair. 

For J/t = 0.4 and 1/8 doping, the /c-dependence of the zero frequency gap 4>(k x , k y , to — 0) 
is plotted in Fig. 1 for the bonding (k y = 0) and antibonding (k y = tt) bands. For comparison, 
the solid curves correspond to (fi (a cos k x — cos k y ) with (fi = 0.3 and a = 0.8. This d-wave- 
like /c-dependence of the gap function is similar to the behavior found in previous studies of 
both the t-J 8] and Hubbard ladders 0]. It implies that the spatial structure of the pairing 
interaction is dominantly a near-neighbor interaction. 

The frequency dependence of the real and imaginary parts of <p(k, to) are shown in Fig. 2 
for k values which are near the fermi surface of the bonding (red) and antibonding (green) 
bands respectively. In general, the gap function 0(fc, uj) is a complex frequency-dependent 
function <p\ + i<p2 with the imaginary part associated with dynamic decay processes. For our 
doped 2-leg ladder with J/t = 0.4, the magnitudes of the zero frequency gaps A = A(kp, 0) 
obtained from eq. (|5|). for both the bonding and antibonding fermi points, are of order 0.15t. 
The gap, A , is reduced from 4>(k F , 0) due to the renormalization factor Z. The onset of 



the imaginary part of <p2{k,uj) seen in Fig. 2(b) appears to occur somewhat below 3Ao. In 
addition, the peaking in (j) 2 (k,uj) and the rapid rise in <pi(k,uj) as this onset frequency is 
approached suggest that a particular excitation mode occurs at a frequency Q such that 
uj = A + fi determines the onset of faik, 

The occurrence of an onset peak in (j) 2 (k,u) at threshold and the short-range nature of 
the interaction involving scattering of pairs from the bonding to antibonding band (q y ~ 7r), 
imply that the pairing in this energy regime is mediated by an S = 1 channel. This follows 
from the form of the coherence factor which varies as ~ ( 1 ± ^%r^r4i5r-^ I , with the plus 

2 y E(p+q)E{p) J ' V 

sign associated with the charge channel and the minus sign with the S = 1 spin channel. For 
a "d- wave-like" gap with q = (hp (bonding) — kp (antibonding), tt) and uj near threshold, 
the coherence factor goes to 1 for the spin channel and vanishes for the charge channel. 
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solid diamonds show a collective S = 1 bound magnon-hole-pair [15J, 3] mode and the open 
symbols the S = 1 particle-hole continuum. A measure of the spectral weight of the 5* = 1 
channel which couples to (fi(k,u) is given by the d-wave projection of the spin fluctuation 
spectral weight 

V d (uj) = — y^(cos k x - cos k y ) (cos k' x - cos k' v ) S{k - k', uj) (7) 

k,k' 



with 

SM^-llmfa — ^1— - S-). (8) 

Here 5° is the Fourier transform of the a spin component. A plot of Vd{oS) versus uj is shown 
in Fig. 3(b). The peak at Q ~ 0.15t arises from the bound magnon-hole-pair mode and the 
high frequency weight comes from the particle-hole continuum. We believe that the onset of 
(p2{k,uj) seen in Fig. 2 is due to a process in which a single particle excitation with energy 
uj = Q m + A emits a bound magnon-hole-pair and drops down to the gap edge while the 
remaining (p2(k,uj) ^ region arises from a coupling to the S = 1 continuum. 

In order to further characterize the dynamic nature of the pairing interaction, we intro- 
duce 

Ift f!)=lf^«M. (9) 

Jo UJ' 

Then, since the gap function satisfies a dispersion relation, one has 



(k, 0) = I(k, n -> oo) + st atic(A;) . (10) 



Here, s tatic(^) represents a non-retarded contribution. For example, if one were to make a 
mean-field approximation in which 



one would have an effective attractive pairing interaction which is independent of frequency 
which would contribute to Static- In Fig. 4(a), we have plotted I(k,Q)/<fii(k,0) versus Q 
and one sees that for J/t = 0.4, I(k, Q) saturates at over 80% of 4>i(k, 0). This implies that 
the dominant part of the pairing interaction comes from dynamic processes. Also shown, 
Fig. 4(b), are similar results for J/t=0.8. In this, unphysically large, J/t regime only 60% 
of the u = gap function is associated with a dynamic pairing interaction. This is similar 
to earlier results for the spatial structure of a pair Q] in which near-neighbor sites for 
the two holes making up a pair were favored for large values of J/t. For physical values of 
J/t, next-near-neighbor (diagonal) hole-hole occupation was favored. The diagonal structure 
arises dynamically and reflects the retarded nature of the pairing interaction for physical 
values of J/t. 

In summary, knowledge of the k- and tu-dependence of the gap function provides informa- 
tion on the spatial and temporal structure of the pairing interaction. Numerical solutions 
for the two-leg ladder show that the pairing interaction has a short-range, near-neighbor 
form, leading to momentum transfer processes which scatter pairs between the bonding and 
antibonding states. The d-wave-like momentum dependence of the gap function and the uj 
onset of 4>2{k, (jS) imply that the pairing interaction in this energy regime is mediated by the 
5 = 1 channel. This channel contains both a bound magnon-pair state and a continuum of 
particle-hole excitations. We believe that the magnon-pair mode is responsible for the onset 
behavior seen in (f>2(k, u). The dispersion relation for </>i(fc, 0) shows that for physical values 
of J/t, the dominant part of the interaction is dynamic with contributions coming from both 
the magnon-hole-pair mode and the particle-hole continuum. 
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FIG. 2: Real (a) and imaginary (b) parts of the gap function 4>(k,cu) versus u for a 2 x 12 ladder at 
values of k near the bonding (red) and antibonding (green) fermi momenta [~ (-^, 0) and ~ (f , tt) 
respectively] for J/t = 0.4. Here, Ao ~ A(kp,AQ) and tt m is the minimum energy for a bound 
magnon- hole-pair excitation. 
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FIG. 3: (a) Low-energy electron- hole S = 1 excitations for a 1/8-doped J/t = 0.4 ladder for tj^ = 
(circles) and q y = ir (diamonds) momenta. The solid diamonds denote a bound magnon-hole-pair 



collective mode 



(b) "d-wave" projection of the spin fluctuation spectral weight with Q m the 



minimum energy of the bound magnon-hole-pair collective mode (1/6 doping). 
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FIG. 4: The integrated contribution of (j)2(k,uj)/uj to (pi(k,0) for the bonding (red), antibonding 
(green) fc-values shown in Fig. 2. (a) for J/t = 0.4 and (b) for J/t = 0.8. The non-retarded 
contributions are shown by arrows. 
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